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P^ Abstract. The magnetic field of the quiet-Sun chromosphere remains a mystery 

^^ for solar physicists. The reduced number of chromospheric lines are intrinsically hard 

to model and only a few of them are magnetically sensitive. In this work, we use 

^ a 3D numerical simulation of the outer layers of the solar atmosphere, to asses the 

I reliability of non-LTE inversions, in this case applied to the Ca ii /18542 A line. We show 

O that NLTE inversions provide realistic estimates of physical quantities from synthetic 

^ observations. 
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jy^ During the past decades, the Ca ii lines have been used to study the solar atmosphere. 

-^ The strong H & K lines and the infrared triplet have been particularly popular in 

• chromospheric studies (see Shine & Linsky||1974t Rutten & Uitenbroek||1991 , Socas- 

O '^' ~-- 

(N 



3 00at 
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2009^ Rutten_20irj and references therein). 



Navarro et al."20 00at [Lopez Ariste et al.||2001HPietarila et al.||2006HLeenaarts et al. 



The interest for the infrared lines has grown during the past decade as un-precedented 
observations wqtq carried-out (Cauzz i et al.||2008| ) and new diagnosing tools were de- 



veloped ( [Socas-Navarro et al.„2000b^ . During the past 10 years, the development of 

?N a new generation of Fabry-Perot interferometers has contributed to the acquirement of 

C^ high-resolution observations over a relatively large field-of-view on the Sun, with a 

restricted wavelength coverage, using the Ca ii /18542 A line, making chromospheric 

studies even more attractive. 

To achieve high polarimetric sensitivity, long exposure times are needed. How- 
ever, chromospheric dynamics limit the the exposure time in chromospheric observa- 
tions, especially if image reconstruction is needed ( |van Noort & Rouppe van der Voort 



2006| ). To illustrate this problem, a time series acquired with 5 seconds cadence is il- 



lustrated in Fig. [T| The difference between consecutive images is large enough to show 
chromospheric motions in very small time scales. 

Additionally, the relatively low effective Lande factor geff = 1 . 1 of the Ca ii /18542 A 
line, and the wide Gaussian core (a common feature in all the lines of the Ca ii IR triplet) 
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Figure 1. Time evolution of a 10.3'' x 7.2'' patch of quiet Sun at the core of the 
Ca II /18542 A Hne. 



gives low polarization levels compared to the narrower Fe i lines, commonly used in 
photospheric spectropolarimetry (geff > 2). 

Unfortunately, scattering polarization also must be considered to interpret quiet- 
sun observations ( [Manso Sainz & Trujillo Bueno|2010HTrujillo Bueno,2010^ . For these 
reasons, most attempts to measure, reconstruct or understand the magnetic topology of 
the chromosphere have focused on active regions and solar network (|Socas-Navarro 



et al. 1, 2000a, ,L6pez Ariste et al. 2001, Pietarila et al. |2007b) and off-limb spicules 



( [Trujiilo Bueno et aL|2005^ [Lopez Ariste& Casini|2005[|Centeno et aL|20T0| ). 

The Ca ii /18542 A line cannot be modeled assuming local thermodynamic equi- 
librium (LTE hereafter). However, the statistical equilibrium equations (one of the sim- 
plest non-LTE recipes) can be used to compute the model-atom population densities 
( |Wedemeyer-Bohm & Carlsson|201 1| ). Furthermore, in active regions, it is reasonably 
safe to assume Zeeman induced polarization. [ Socas-Navarro et aL]pOOOb[ ) presented a 
scheme to carry-out NLTE inversions of spectropolarimetric data, assuming statistical 
equilibrium and Zeeman polarization. 



Two datasets acquired with the Swedish 1-m Solar Telescope (SST, |Scharmer et al. 
2003| ) and the CRIsp SpectroPolarimeter (CRISP, Scharmer|2006| ), are described below 
illustrating the complexity of the observational data. 

Figure [2] shows very quiet-Sun data at -231 mA from the core of Ca 11 /18542. The 
clapotisphere described in |Rutten| ( |2011| ) is clearly visible all over the field-of-view. 
Fibril-like features originate from bright points located in the center of the image. The 
weak signal in the Stokes V image originates in bright points and large-scale patterns 
are not visible. The signal in Stokes Q and U is, as expected, below the noise level. 

Active regions show a completely different scenario. The Stokes Q, U and V pan- 
els in Fig. |3] show areas with strong signal. The complicated features seen in Stokes V 
are not all due to the magnetic field topology, but to Doppler motions and emission in 
the Stokes / profile, which can lead to a reversal of the V profile. 

The work presented here attempts to asses the reliability of NLTE inversions and 
how to interpret the inversion results. For that purpose we have created synthetic ob- 
servations from a 3D simulation of the solar chromosphere, evaluating the radiation 
field in 3D. Those synthetic profiles are degraded with Gaussian instrumental profiles 
and random noise before inverting them. The resulting magnitudes are compared with 
those from the original 3D snapshot. 

NLTE inversions have been used in the past to infer physical quantities in the solar 
chromosphere (see [Socas-Navarro et aL]|2000b[ [Pietarila et al.[[2007a| , however, the 
work partially presented in this paper aims to asses the reliability of the results and 
evaluate possible systematic errors induced by the physical assumptions used during 
the inversion process, some of them listed below: 
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Figure 2. Very quiet Sun observation at -23 1 mA from the core of the Ca ii /18542 
line. Stokes / and V are shown on the left and right panels respectively. The data 
were acquired on 2009-05-29 and the field size is 10'' x 10'\ 



• 3D NLTE efFects: the radiation field propagates through the 3D structure of the 
solar atmosphere, however, the inversion decouples it, treating each column as a 
plane-parallel atmosphere (usually referred to as 1.5D). 

• Hydrostatical equilibrium and the ideal gas equation of state are assumed to al- 
low physical consistency between electron pressure, electron density and temper- 
ature. 

To evaluate these two eflFects, we create realistic synthetic observations from a 
snapshot from a 3D numerical simulation from the photosphere to the corona. The in- 
versions are carried out using our synthetic observations, previously degraded assuming 
a list of instrumental profiles and noise. 

2. Radiative transfer 



This work is based on that carried out by [Leenaarts et al.| ( |2009| ). To create realis- 
tic observations we use a 3D simulation computed with the Oslo Stagger Code (see 
Hansteen et al.|[2007| ). The snapshot extends from the upper convection zone to the 
corona, and covers a physical patch on the surface of the Sun of 16.6 x 8.3 x 15.5 Mm 
(256 X 128 X 160 pixels). 

The snapshot has been interpolated to the higher depth-resolution but restricted to 
the zone of importance for the radiative transfer calculations. The simulation has an 
average magnetic field strength of 150 G in the photosphere. 

The NLTE radiative transfer calculations are performed with a Ca ii model atom 
consisting of 5 bound levels plus a continuum. The synthetic observations are com- 
puted in two steps. To compute the populations densities we use the radiative trans- 
fer code MuLTi3D ( [Leenaarts & Carlsson||2009| ) that allows to evaluate the 3D radia- 
tion field (3D NLTE hereafter). The resulting population densities are used by Nicole 
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Figure 3. Active region observed at -225 mA from the core of the Ca ii AS542 
Hne. Stokes /, Q, U, V are shown from left to right and top to bottom respectively. 
A photospheric image is shown in the central panel to set the context, where the tick 
mark separation is W and the total FOV 53'' x 53''. Data processing by J. de la 
Cruz Rodriguez and L. Rouppe van der Voort. Observations by A. Ortiz-Carbonell 
andB. Hole (2011-09-24). 
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Figure 4. Synthetic observations at 209 mA from the core of the Ca ii /18542 A 
Hne. The Stokes Q and U images are scaled between ±6 • 10"^ and Stokes V is scaled 
between ±1.5 • 10"^, relative to the continuum intensity. The white box indicates the 
area where inversions are carried out. 



(see H. Socas-Navarro contribution) to compute the full-Stokes vector. The full-Stokes 
profiles are normalized to the average continuum intensity. 

Figure |4] shows the resulting Stokes /, Q, U and V images at -209 mA from the 
core of the Ca ii /18542 A line. 

The synthetic spectral profiles are convolved with Gaussian instrumental profiles 
of full- width-half-maximum FWHM = 50 mA. Additive photon noise is added after- 
wards with standard deviations cr = [10""^, 10"^]. 

In Fig. |5} the Stokes Q image from the central part of our snapshot is shown at 
-1-205 mA from the line core. The profiles have been convolved with Gaussian instru- 
mental profiles of full- width-half-maximum (FWHM hereafter) AA = 50 mA. Random 
photon noise has been added with standard deviation cr = [0, 10""^, 10"^]. Most spectral 
features in Stokes Q and U have signals of the order of 10""^, peaking at most at 10"^ 
where the magnetic field is strong. 

In the middle row, corresponding to cr = 10""^, most polarization features are 
above the noise level, regardless of the spectral resolution. However, in the bottom row 
(cr = 10"^), the limited spectral resolution and noise wash all the features from the field 
of view. 

From an observational perspective, the magnetic field orientation and inclination 
are computed using Stokes Q and U. To infer the full magnetic field vector, it is crucial 
to carry out observations with noise well below 10"^, as it is shown in Sect.p] 
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Figure 5. Synthetic Stokes Q images computed at +205 mA from the core of 
the Ca II /18542 A line. The spectra have been degraded with Gaussian profiles 
and random photon noise. From left to right, the full- width-half-maximum of the 
Gaussian profile increases from AA = mA (no convolution) to AA = 100 mA. 
From top to bottom, the noise is increased from cr = to cr = 10"^. The images are 
scaled to ±0.15%, relative to the continuum intensity. 
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3. Inversions 

The inversions are carried out using the inversion code Nicole. The results shown in 
this section correspond to chromospheric averages computed from Log(T5oo) = -4 to 

Log(T5oo) = -6. ^__^ 

The inversions are initialized using the VAL-C model ( Vemazza et aL|198T] ). Dur- 



ing the inversion, the variables of the model are modified at node points placed equidis- 
tantly along the depth scale, to fit the observed profiles. 

Each inversion is computed in two cycles, a similar approach to that described 
by |Pietarila et al.| ( [2007a| ); |Socas-Navarro| ( |201l] ). The first cycle aims to search for a 



rough solution using few node points. In the second cycle more degrees of freedom are 
allowed to refine the solution. The inversions are computed using critically sampled 
spectral profiles, assuming that the spectral resolution is given by the FWHM of the 
instrumental profile. 

The results from the inversions are shown in Fig. [6j The reconstructed B^ is less 
noisy than |B;c| and |B^|, an expected results given that B^ is determined from Stokes V, 
where signals are much stronger than in Stokes Q and U. 

Noise is a critical factor for the inversions. The inversion does a reasonably good 
job reconstructing B^ for all the noise levels considered in this work. The diff'erences ap- 
pear with |B;f| and |B_y|. The results obtained with cr = 10""^ are similar to the noise-free 
case where there is good correspondence between the inversion and the 3D simulation. 
At 10"^, all the magnetic features are lost, under the noise level. 

The inverted temperature and line-of- sight velocity (l.o.s. hereafter) are insensitive 
to the spectral degradation and noise levels considered here, thus only one panel is 
shown in Fig. [7] The l.o.s. velocity is recovered mostly from the Stokes / and V 
profiles. The inversion shows similar feature to those in the 3D model, with minor 
diff'erences probably due to diff'erences in the quality of the fit and the rough interval 
that has been selected to compute the depth averages. 

The inferred temperature map has less contrast than the 3D simulation. We ex- 



pect 3D NLTE eff'ects to work in the same direction. |Leenaarts et al.| ( |2009| ) reports 
decreased contrast in intensity monochromatic images computed in 3D NLTE when 
compared with the 1.5D case. This eff'ect would lead to decreased contrast in temper- 
ature when inverting those profiles, however it is hard to quantify by how much. The 
hydrostatic equilibrium assumption does not seem to have a large impact on the inverted 
quantities. 



4. Conclusions 

In this paper we present preliminary results from an ongoing work that will be presented 
in detail in a forthcoming publication. 

We have created synthetic Ca ii /18542 observations in NLTE, evaluating the 3D 
radiation field and assuming Zeeman induced polarization. The resulting profiles have 
been degraded with a list of instrumental profiles and noise and inverted afterwards. To 
infer the full magnetic field vector from observations, the noise level must be of the 
order of 10""^ and clearly below 10"^. 

The inversions are computed assuming plane-parallel geometry and imposing hy- 
drostatic equilibrium. Despite these assumptions, the NLTE inversion is able to sue- 
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Figure 6. Comparison of the magnetic field vector from the 3D snapshot and the 
results from the inversion. Gaussian instrumental profiles with FWHM = 50 mA 
have been used. From left to right, the panels show B2(x, j), |Bx(x,3;)|, |By(x, j)|. 
From top to bottom, the magnetic field from the 3D snapshot, followed by the inver- 
sion results assuming random photon noise of cr = 0, 10""^, 10"^ correspondingly. 
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Figure 7. The results from the NLTE inversion are compared with the quantities 
present in the 3D simulation. The top row corresponds to temperature and line-of- 
sight velocity from the 3D model, and the bottom row to the inversion results. These 
panels correspond to a chromospheric average from Log(T5oo) = -4 to Log(T5oo) = 
-6. 



cessfully reconstruct the chromospheric quantities from the 3D simulation, although 
the quality of IB^^I and |B_y| strongly depends on noise. 

The 3D NLTE effects mentioned in Sec.|3j tend to decrease the contrast of intensity 
profiles. The decreased contrast in the inverted temperature map could be explained 
with this eff'ect. 
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